Universality class of the restricted solid-on-solid model with hopping 
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■ We study the restricted solid-on-solid (RSOS) model with finite hopping distance Iq, using both 
analytical and numerical methods. Analytically, we use the hard-core bosonic field theory developed 

£S) . by the authors [Phys. Rev. E 62, 7642 (2000)] and derive the Villain-Lai-Das Sarma (VLD) equation 

- .. ' for the Iq = oo case which corresponds to the conserved RSOS (CRSOS) model and the Kardar- 

Parisi-Zhang (KPZ) equation for all finite values of Iq. Consequently, we find that the CRSOS model 
belongs to the VLD universality class and the RSOS models with any finite hopping distance belong 
to the KPZ universality class. There is no phase transition at a certain finite hopping distance 
contrary to the previous result. We confirm the analytic results using the Monte Carlo simulations 
for several values of the finite hopping distance. 

'"rH 1 | PACS number(s): 05.90.+m, 81.10.Aj 
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1 ■ I. INTRODUCTION 

In recent years, the field of nonequilibrium surface growth has been investigated using various discrete models and 
, continuous equations Q . The comprehension of nonequilibrium surface growth plays an important role in understand- 
ing and controlling many interesting interface processes, such as vapor deposition ||, crystal growth molecular 
beam epitaxy (MBE) [Q, and so on. During the MBE growth process, the conserved growth condition is applied 
without defects, such as overhangs and vacancies. Various discrete conserved models for MBE, describing the kinetic 
properties of this type of surface growth, have been proposed and studied by intensive numerical simulations. 

The main purpose of studying discrete models is to measure scaling exponents for the kinetic roughening, which 
£j ' determines the asymptotic behavior of surface growth on a large length scale in a long time limit. The important 
result of the kinetic roughening studies is that a large variety of different discrete growth models can be divided into 
only a few universality classes. The surface width W, which measures the root-mean-square fluctuation of the surface 
height, scales as 

C*' W(L,t)~L a f(t/L*), (1) 

(N : 

£NJ . where the asymptotic behavior of the scaling function f(x) is constant for x ^> 1 and x p for x <C 1 with /3 = a/z. The 
scaling behavior of the growth is characterized by three exponents: the roughness exponent a, the growth exponent 
C^j , /3, and the dynamical exponent z. These exponents determine the universality class. 

In the coarse-grained picture, evolution of the growing surface is usually described by a stochastic differential 
equation (SDE) for the height variable /i(x, t) as a function of the surface coordinate and time. For discrete models 
of MBE growth, several SDEs were suggested and it was generally believed that there is a correspondence between 
discrete growth models and continuum SDEs. The common way of establishing the link between discrete models 

■ and continuous equations is a simple comparison of critical exponents determined from computer simulations of the 
' discrete model with exponents for the continuous equation. 

There have also been attempts to establish the correspondence in an explicit way. The systematic procedure for 
O ■ establishing a continuous equation corresponding to discrete models, starting from the master equation in discrete 
J> \ space was proposed by Vvedensky et al. || and has been successfully applied to the derivation of growth equations for 
T '"j ■ some discrete models, including the solid-on-solid (SOS) model, the restricted solid-on-solid (RSOS) model, as well as 
the Wolf- Villain and Das Sarma- Tamborenea models |^-||]. However, there are several difficulties with this procedure; 
in particular, in converting from the equation system for a discrete set of heights to an equation for a continuous 
function /i(x, t), the procedure requires the regularization step, in which the nonanalytic quantities are expanded and 
replaced with analytic quantities, i.e., the step function is approximated by an analytic shifted hyperbolic tangent 
function expanded in a Taylor series. However, the form of the regularized function is uncertain, and different choices 
of this function lead to different results. Some forms of the regularized function have been suggested, but the problem 
of a proper choice of a regularization scheme has not been discussed. 

To overcome this kind of uncertainty, the authors proposed a new schematic formalism M , deriving the continuous 
equations, such as the Edwards- Wlikinson (EW) jl0| and Kardar-Parisi-Zhang (KPZ) [pd| equations, from the body 
centered solid-on-solid model and the RSOS model. In this paper, we apply our formalism to a new kind of MBE 
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growth model proposed by Kim, Park, and Kim |l2| . This model allows the deposited particle to relax to the nearest 
site where the RSOS condition on neighboring heights is satisfied and has the conserved growth condition constraint, 
which means the deposited particles are possible to hop for an infinite distance until they eventually find the site with 
the RSOS condition satisfied. Applying our formalism to the above conserved RSOS (CRSOS) model, we not only 
derive the Villain-Lai-Das Sarma (VLD) equation |13] for the model which belongs to a different universality class 
from the EW and KPZ equations, but also we are able to predict the coefficients in the VLD equation, which was not 
possible by other methods. 

Observing that the RSOS model belongs to the KPZ class and the CRSOS model belongs to the VLD class, we 
went one step further to study the RSOS model with the finite range hopping (RSOS/H). In this RSOS/H model, 
it is possible for the deposited particles to hop a finite distance Zo until they find the site with the RSOS condition 
satisfied. If they fail to find the site with the RSOS condition satisfied within the distance Iq in both directions, the 
deposition process is rejected. The RSOS model corresponds to the RSOS/H model with Iq — and to the CRSOS 
model with Zo = oo. We apply our formalism to the RSOS/H model with Zo finite and find that this model belongs to 
the KPZ class, contrary to the previous report by Kim and Yook fLU , who concluded that there is a phase transition 
between the KPZ class to the VLD class along the parameter Zo. 

In Sec. H our formalism to derive the continuous equation from the discrete model is briefly explained and the 
procedure of derivation is described. The detailed calculations are attached in the appendices. Extensive numerical 



simulations are presented in Sec. Ill and the summary and discussion are given in Sec. [V 



II. DERIVATION OF THE STOCHASTIC EQUATION 

In this section, we derive continuous equations for the one-dimensional RSOS/H model with a hopping distance Zo 
and for the CRSOS model corresponding to Zo = oo. We restrict ourselves to the case wherein the height difference 
between two nearest neighbors is not larger than 1. For a succinct description of the dynamics, we introduce the 
nomenclature that if a site m satisfies the condition \h m + 1 — h m ±i\ < 1, this site is called "stable." Following this 
nomenclature, the growth algorithm of the one-dimensional RSOS/H model is as follows: (i) We choose a site m 
randomly, (ii) Sites from to — Zo to to + Zo are examined to determine if they are stable sites, (iii) If a stable site is 
found within the interval from to — Iq to m + Zo, a new particle is deposited to the nearest stable site from to (to itself 
can be a candidate for deposition). However, if stable sites are nonexistent in the examined interval, the particle drop 
is rejected and the system remains unchanged. After this try, the time is increased by l/L, where L is the system 
size. We assume periodic boundary conditions. 

Since the height difference between two nearest neighbors is restricted not to be larger than 1, the RSOS/H model 
is mapped onto the reaction-diffusion system of hard-core particles with two species. The step-up (-down) is mapped 
to an A (B) particle. If two nearest neighbor sites have equal height, a particle vacuum is located between these 
two sites. The site where the particles reside is labeled by an integer, and the site for height by a half-integer. This 
mapping is depicted in Fig.|l[ According to this mapping, the dynamics of the RSOS/H model can be described by 
the (imaginary time) Schordinger equation (d/dt)\^f; t) = —H\^; t) for the state vector \9\ t) = J2 C P(C; t)\C), where 
P(C;t) is the probability with which the system is in state C at time t, and H, called a Hamiltonian, is an evolution 
operator 

H = (f ~ R n - An) Ln, (2) 

n 

where 



j=l \k=l fc=l / 



(3) 



and 

R n = a n a n + b n+l b n+ \ - a n Z>Jj +1 & n+ i, 

A n = [a n + b n J (a„+i + bl +1 ) . 

The role of the rejection operator R n is to check the stability of site rt + i , that is, if a configuration \C) has a stable 
site at n + i, R n \C) — and otherwise R n \C) = \C). The adsorption operator A n describes the configuration change 
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after a successful deposition. a n (b n ) is the annihilation operator of an A (B) particle at site n and a n (b' n ) is the 
corresponding creation operator, satisfying the mixed commutation relations presented in Ref. |J, 



Q'nQ'ri — bnb n 



I 



a n a n 



btJ> n , 



ilbl = 0. 



(5) 



All operators at different sites commute with each other. 

To find the SDE for the RSOS/H model, we apply the method recently introduced by the authors ||. First we 
assume the existence of the lattice version of the SDEs in terms of density. Those equations arc supposed to take the 
forms 



dt a n = Ct ({a}, {&}) + J2 bn^iit) + gi B n C{t)} , 

m 

d t b n = C* ({a}, {&}) + £ [g^C(t) + 9X(t)] , 

m 

where dt = d/dt and , ^ are white noises with correlation 

(X,X' = cither A or B), 



(6) 
(7) 



(8) 



where S nm and 6x,x' are Kronecker deltas and S(t — t') is the Dirac delta function. The matrix g is related to the 
Kramers-Moyal coefficient in such a way that 



Er,X»9ir X 9*r X = C™ ({a}, {b}) 

(X, X', X" = either A or B). 



(9) 



Here we are using the Ito interpretation. The field a (b) in the curly bracket represents the density of species A (B) 
at all sites. From here on, without a hat above itself a mathematical symbol is a pure number as opposed to an 
operator. a n should not be confused with a n . The former is a density at site n that runs over real numbers, while 
the latter is an annihilation operator. By requiring that the noise average of observables in Eqs. (^|) and ([?]) has the 
same behavior with the ensemble average of the number operator, we find Q 



(10) 



<c»>* 


= ([H,aia n }) t , 






= {{H,b\b n ]) u 




(fiAA\ 


= ([aia n , [H, a] n a m ]]) t , 




— (@mn)t — ^ 


ala n , [H, b\J) m } 




= (jbib n ,[H,hib m ]]) t , 



where the (• • -)t on the left hand side represents the average over noise at time t and that on the right hand side 
stands for the ensemble average. The arguments of the Kramers-Moyal coefficients are dropped for brevity. 

As presented in Ref. |)| , the ensemble average of any operator can be interpreted as an average of number operators 
due to the property of the projection state (-|, which is defined as a sum over all possible microscopic configurations 
and is itself a left eigenstate of H with eigenvalue 0, and in turn, the ensemble average of a number operator is 
mapped to the noise average of density. This procedure leads us to find the Kramers-Moyal coefficients C„ , , in 
terms of the density fields. We call this procedure figurization, which means "expression in number." To represent 
the figurization, we use the symbol i — > to the left of which is an operator (or a product of operators) and to the right 
of which is the corresponding density representation. 

To complete the derivation, we perform the commutation relations between the Hamiltonian and the density 
operators and so forth: 



H, a n a n 
alai, H,d)a n 



= a} n (a n+ i + b* n+1 j L n - a n [a'l_ 1 + 6„_i j t n -\ 
= -b n (a n+ i + b n+1 ^j L n + bl f aj,_x + K-i^j L n 
= Sin ( al ( a n+ i + P n+l ) L n + a n ( a 1 n _ 1 + 6„_i ) L n 



(11) 
(12) 
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ajai, 



H, b n b n = Sin (b n (a n+ i + b\ + ^j L n + P n (a n _ 1 + 6 n _iJ £„~i) 
H, b n b n = 5i^ n+ ib n a n+ iL n + 5i )n _iSj [ a n _ 1 -L ra _i, 



(13) 

(14) 
(15) 



where S n i is a Kronecker delta. Following the figurization, we find the Kramers- Moyal coefficients. For a later purpose, 
we give some examples of the figurization. The figurization of R n is R n i — > R n = a n + 6 n +i — a n b n+ i , and the symbolic 
representation of the figurization of the product of i?'s is 



fe=i 



n-\-k 



R 3 



(16) 



where the superscript j should not be confused with the power. When j = 1, i?^ is denoted as R n +i- R 3 n are given 
by the following recursion relations: 



j'+i 



R 3 



R 3 n = a n+1 R 3 n+ \ + (1 - a n+ i) ]^[ b n+k , 

k=2 

i+i 

&n-l-R^Ij_2 + (1 - bn-l) ]^[ a„_fc, 



(17) 
(18) 



fc=2 



where j > 1 and we define = 1. The physical meaning of Eqs. (|T^) and ( p8| ) is as follows: We divide the situation 
that prohibits the height increase at sites n+|,n+|,---,n + j -f-T by the condition at site n + 1. If there is an A 
particle at site n + 1, the height increase is suppressed at site n + | irrespective of the condition at site n + 2. Hence 
the first term of Eq. ( |l7j ) follows. If there is no ^4 particle at site n + 1, to suppress the deposition at site n + §, 
there must exist a _B particle at site n + 2 and this should be continued until the site n+j + 1, because at site n + k 
(2 < k < j) no A particle is present; this condition is represented by the second term. To comprehend Eq. (p"8|), we 
only have to perform the mirror transformation relative to site n. By the mirror transformation relative to a site n, 
we mean the exchange of a and b (<z <-> b), followed by n + k <-> n — k. Under this transformation, R 3 n+k changes into 

R 3 n _ k _j_2 for an arbitrary k. The mirror transformation of Eq. ( p"7j ) is Eq. (p"§|). 

With these notions, we will find the SDE of the RSOS/H model. At first, the deterministic part of SDE is found. 
Since the main concern is not the respective dynamics of the A and B particles, but it is D n = a n — b n (the local 
slope) and S n = a n + b n (the slope density), we will write the SDE for D and S rather than for the A, B particles. 
The Kramers-Moyal coefficient is obtained by subtracting Eq. ( p"2] ) from Eq. (11) and C„ by adding Eq. ( p"2] ) 
and (0) followed by the figurization. 



CB = 



1 2l ° r 



K-2 - +K~ (RLj-2 



1R 3 



R 3 



2io + l 



2 ( R n-2 



p2i + l 



R 



2l + l 



R 



2/„ 



n-1 ' ±l -n-2l -2 1L n-2l -lJ ' 

3(o„ + b n ) — {b n a n+ i + a n b n ^i — 36„6„+i — 3a„a„_i) 



(19) 



1 / \ 

- (K - K-2 - K-j + K-j-2) - bn(l - S n+1 )b n+2 - On(l - 5 n _x)a„_ 

l/p2/o + l , p2io+l , R 2; + l i R 2; + l N 
Q^ri-l ' /x n-2i -l ^ «-2 ' ^71-210-2) 

j+2 

b n a n+ iR 3 n+1 + b n (l - S n+ i) K+k + M 1 - K+l)R J n -j- 

k=2 
j+2 

a„6„_ii?^ l _ :( -_ 3 + a„(l - S n -%) Y[ a n-k + a n (l - a n -i)R 3 n 



2 

2(0-1 

E 

J=i L 

2io-l 

E 



fc=2 



(20) 
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As pointed out in Ref. [ JL5| , the mass term in Eq. (|2(j) makes the step density saturate fast. As a result, the step 
density S becomes a slave field of the slope D and takes the form 



S = 2p(l ) + n (l Q )dD + 9 (Iq)D 2 + 



(21) 



where 2p(Zo) is the (mean- field) stationary step density and pq and 9q may depend on Iq. Here d = (d/dx). Since 
this system has the mirror symmetry whose continuum version is invariant under the transformations D — ► —D and 
d — > —9, wc do not expect the occurrence of D in Eq. (|2l"|). The parameters p(Zo), Mo('o)> an d 0o('o) are determined 
by the stationarity of Eq. (pp|). In Appendix |A|, we obtain p, 9q, and as functions of Iq and show the approximate 
solutions of these parameters. We can now rewrite Eq. ( |l9| ) in terms of the height field h using Eq. (^) and 



D = dh. 



(22) 



The last task is to determine the noise strength. This is accomplished using Eqs. (||) and (R). After eliminating S in 
favor of D, there is only one kind of noise with strength 



c — n 



AA 



^iBA 



Using figurization of Eq. ( |10[) and keeping the most relevant term, we find 

C nm — —(1 — p) 2 L(5 n ^ m +i — 2S nm + 8n t m-l) + 



where 



For finite Iq wc found 



L = J^p l (l + l-pl). 



1=0 



(23) 



(24) 



(25) 



where 



d t h = Voc + vd 2 h + -(Oh) 2 + £(x, i), 
(Z(x,t)ax',t'))=D ii 5(x-x')5(t-t'), 



v ao = l- R{2l a + 1) ~ 1 - p\ 



oo 



v = -p(2l + 1) ~ p»°(l + 1)(2Z + 1) + + O (p£>) 



A = -20(2/ o + 1) ^ + Io)(l + 2Z ) 1 + 



6 4 
V2 + 1 



h +O(p« ) 



(26) 



(27) 



2V2- 1 



1-P 



,. , 6 + 5V2 19V2-22, 
1 1 Z 



2<S 
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with ^ = (2 - \/2)/2. The numerical values of p(lo), v(Iq), and A(Zo) are given for several values of Iq in Table I. 
For infinite Iq 



where 



dth = woo - vd A h + \d 2 {dh) 2 + £(x, t), 
(i(x, t)i(x', t')) = D (( 6(x - x')S(t - t'), 



Voo = 1, 

. 21-12\/2 



(28) 



A 



10-3\/2 
2 ' 
2V2-1 



(29) 



5 



Equation ( |26| ) directly shows that the RSOS/H model, for any finite lo, belongs to the KPZ class, and Eq. ( |28| ) 
suggests that the CRSOS model is described by the VLD equation. However, the first line of Eq. (|l9| ) has the form 



- 2Ri_ x + Ri - {mirror terms of «_ 2 - 2W n _ x + R 3 n )}}. (30) 

which has the following implications: If the continuum version of -R^_ 2 — ^R-n-i + -R^ has a nonvanishing coefficient of 
d 2 (D 2r+1 ) with a non- negative integer r, it is not certain that Eq. (|30| ) is a lattice version of <9 3 ; consider the mirror 
transformation in the continuum limit. The occurrence of d 2 (D 2r+1 ) in Eq. ( |30| ) is directly related to the appearance 
of d(dh) 2r+1 in Eq. (|2§|), which generates an EW term by the dynamic renormalization group fl7fl . Appendix [b| 
shows that this is not the case. The vanishing of D 2r+1 in the continuum limit of R J n+k guarantees the vanishing of 

d 2 (D 2r+1 ) in R 3 n _ 2 — 2i?^_ 1 + R 3 n . In view of this, we conclude that the continuum equation of the CRSOS model is 
the VLD equation. 



To confirm this conclusion, we have performed Monte Carlo simulations as outlined in Sec. Ill 



III. NUMERICAL STUDY 



Although the RSOS/H model was studied numerically by Kim and Yook pjj , their results are contradictory to our 
derivation. As a result, we need to perform extensive numerical simulations to verify our results. In the derivation, 
we found that the coefficients of the EW terms and the KPZ terms are vanishingly small, though finite, for large 
lo- Hence, we may find a crossover of roughness exponents from VLD (a v id — 0.95) to KPZ (ak pz = 0.5). We are 
preoccupied with the numerical observation of this crossover. 

In Fig. ^, we draw the saturated width W sa t as function of the system size L for some IqS. The system sizes are 
64, 90, 128, 180, 256, 360, 512, 720, and 1024. For relatively small system sizes, the roughness exponents are near to 
the values reported in Ref. As expected from our derivation, we see a crossover for large system sizes. 

To clarify the crossover behavior, a scaling plot is given in Fig. ||. The anticipated scaling form of the saturated 
width is 



W sat (l ,L)=L a ^g(q/L), (31) 

where 7 is the crossover exponent. The asymptotic behavior of the scaling function g is expected to be 

J const when x — > 00 

ffW ~ j x a vld -a kpz when x _> q • W 

The best fit for the data set shown in Fig. || corresponds to avid = 0.9 and 7 = 2.0. The fitting parameter a v id 
obtained is somewhat smaller than the known value of the roughness exponent of the CRSOS model. This is most 
likely due to the smallness of lo. For example, the data for lo = 10 in Fig. [| yield 0.86 rather than 0.95. The meaning 
of 7 is as follows: When the system size is sufficiently large, we expect the width of the characteristic mountain to be 
~ L Q k PZ D ue t 1^ RSOS condition, the linear size of the mountain is also expected to be ~ L" 1 ^. The smaller the 
system size, the less the rejection events occur due to the shrinking of the characteristic mountain. If lo is comparable 
with the linear size of the characteristic mountain (L Qkpz ), the system starts to behave differently. Eventually, the 
system with small size, L ak ^ <C Iq, cannot feel the existence of Iq- Thus the crossover length L* is expected to be 
?o /akp % that is, 7 = l/a kpz = 2. 



IV. SUMMARY AND DISCUSSION 



We studied the RSOS/H model using both analytical and numerical methods. We derived the continuum equation 
for the microscopic discrete model analytically, and found coefficients of the EW term d 2 h, the KPZ term (dh) 2 , as 
well as the VLD term d 2 {dh) 2 . We observed that the coefficients of the EW and KPZ terms behave as - a 2 lo ll for 
sufficiently large Iq, which is consistent with the previous numerical study Jl4] |. Accordingly, we concluded that the 
RSOS/H model for any finite Iq eventually belongs to the KPZ class and the CRSOS model belongs to the VLD class. 
Numerically, we reported the crossover from the VLD class to the KPZ one, which confirms our derivation. Moreover, 
we found a crossover exponent 7, which is argued to be 2. 
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Besides these studies, we can offer an (nonrigorous piausible) argument to anticipate the universality class of the 
RSOS/H model by employing the block spin concept of Kadanoff 's. Consider a system with linear size L and hopping 
distance Iq. Similar to the block spin in the Ising model, we introduce a coarse-graining parameter b which blocks 
the b sites by one. Although the exact transformation of coarse-graining cannot be determined, we expect that if it 
exists, Iq may be renormalized as ~ lo/b 1 . Hence, we expect that the stable fixed point corresponds to Iq = 0, which 
is the KK model and the unstable fixed point corresponds to Iq — oo, which is the CRSOS model jlq] . 
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APPENDIX A: CONTINUUM LIMIT AND DETERMINATION OF 0o AND p, 

In this appendix, we derive the continuum limit of R l n+k where n is the reference point. One may obtain the 
continuum limit directly from Eq. ([l7]), but we do not follow this route. Rather, we make recursion relation about 
the parameters that appear in front of the field (see below). The use of Eq. (|2l|) makes it possible to represent the 
continuum limit in terms of D : which is directly related to the height by Eq. (|22|). We can find the correct coefficient 
that appears in Eq. (^l|) after we find the continuum limit of R l n+k - 

The continuum limit of R l n+k takes the form 

R l n+k = R{1) + 7 (Z, k)D + 0(1, k)D 2 + fj,(l, k)dD + ■ ■ • , (Al) 

where the implicit dependence on Iq is assumed and the argument of D is dropped for simplicity The symbol "=" 
represents the continuum limit of a quantity. Equation ( |l7j ) allows us to find the recursion relations. Before going 
further, the explicit form of R(l) in the case of no tilt boundary condition is found. The recursion relation becomes 

R(l)=pR(l-l) + (l-p)p\ 
^p- l R(l)=p-V-VR(l-l) + l-p (A2) 

with an (sort of) initial condition R(l) = 2p — p 2 . p is the abbreviation of p(lo) in Eq. (|2l]). It is trivial to find the 
solution that reads 

R(l)= p l + p\l- p)l. (A3) 
To find the recursion relations, we need a continuum limit of , which is defined as 

l+i 

= (1 - a n+k+1 ) Y[ b n+k+m = P^xbn+k+i+i- (A4) 

m=2 

One may directly calculate the continuum limit of Pf , but we follow another path. Let the continuum limit of Pf be 

Pf = (1 - p)p l + 7l (Z, k)D + 0i(l, k)D 2 + m(l, k)8D + ■■■. (A5) 



From Eq. (A4), we find the recursion relations (I > 2) 

7 i(*,fc) _ 7i(i-l,*0 1-P 

p l p 1 - 1 2p ' 

fli(Z.fc) _ 0i(i-l,fc) li(l-l,k) l-p 

p l p 1 - 1 2p l + 2p ' 

Hi(l,k) - 1 - p 

1 — = Pi + {po-k-l- 1 — — , 

p l p l 1 2p 

with initial conditions 



(A6) 
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71(1,*;) = --, 

6i(l,k) = 1(1 + 20 o -Ap0 o ), 



(A7) 



9p a nd /^o are functions of In as shown in Eq. (pl|), which are to be determined by Eq. (p0[). The solutions of Eq. 
(A6) under the condition (A7) are (I > 1) 



7l (/,fc)=^- 1 

fli(z,fe) = P l ~ 1 

f , 1 (l,k)=p 1 - 1 



7l (l,fc)-i_£(l-i; 



1 \ Z- 1 1-p 



(Z - 1)(Z - 2) 



(A8) 



fc) + l-^(/io - fe - 2)(Z - 1) - i-^Z(i - 1) 



Now we can obtain the c ontin uum limit of R l n+k - Using Eqs. (17), (A5), and flAg), we find the recursion relations of 
parameters given in Eq. (Al) 



7 (Z,fe) =p 7 (Z-l,fc+l), 



6(1, k) = p6(l -l,k + l) + 6t(l, k) + yP /_1 [I ~ p(l ~ 1) 

l-i 



with conditions 



p(l, k) = pp(l -l,k + l) + m(l, k) + £— (1 + fc + /i ) [Z - p(l - 1)] 



7 (l,fc) = 0, 
1 



6(l,k) = -+6 (l-p), 



(A9) 



/l(l,fc) = (1 - P ) [HQ - 



1 



We solve Eq. (A9) step by step. It is trivial to find that 7 (Z, k) = 0. As shown in Appendix O, the vanishing of 7 is 
not a coincidence. Since 6%(l, k) has no k dependence, we expect that 6(1, k) also has no k dependence. We find 



6(1) = lp' 



i-i 



'„(] /, ) - ( 20 - 2p0 + — 



1 \ I - 1 1-/9 



24/9 



(/-!)(/ -2) 



(A10) 



where fc is removed from the argument of 6 due to the independence. p(l, k) seems to have an explicit k dependence, 
but by inserting pi(l, k) directly, we find the independence of p(l, k) on k. The result reads (we drop k by the same 
reason as 8) 



M (0 = (i - p)p 



l 2 



(All) 



By requiring C,f = at the stationary state, we can determine p, On, and po as a function of Iq- To find p, Oq, and 
po, we need the continuum limit of 



k=0 



s (l + 4p6 )D 2 + p ll -^-(2p a -l)dD. 



We determine p(ln) from Eq. (|20|). p is the solution of the equation 

2a; 2 - Ax + 1 - x 2 ' 0+1 [l - x(5 + 2l ) + x 2 (3 + 2l )] = 0, 
whose approximate solution is 



(A12) 



(A13) 



8 



p(lo) 



n 2 'n + l/ 



(2/o + l) )+0{ P t°), 



(A14) 



with poo = (2 — \/2)/2, which is the solution for infinite lo- By the same reasoning, we can find 60 and po, whose 
approximate solutions read 



%(lo) - ^ 



1+ ^ofi^! + 21 -^ + 



MM 



4-3V2 



1 + /C - 



4 6 

3\/2 4-17V2 



2-V2 2 %/2 3 

o '0 "I 5~ l o 



12 



, 2 + V2 2 2 + 2V2 3 
'0 H 2 — 'o H 3 '0 



0(p 4 i°) 



(A15) 



For finite lo, the most relevant terms arise from the second line of Eq. (|T^) and the resulting equation is the KPZ 
one, with coefficients 



-p(2l + l), A = -26(2l + 1). 



(A16) 



We give the numerical values of p, v, and A in the table. These numbers are determined from the direct calculation 
of p, 6q and po using Eq. (fl9l). These A's should be compared with the previous numerical results fli)] . Note that A 
in Ref. is one half of A here. 

For infinite lo, the second line of Eq. (|l^) vanishes, of which the physical meaning is that a dropped particle 
eventually find a stable site. Hence the continuum equation becomes the VLD equation with coefficients 



-.00 -.00 

~X>Cj). a = -X>C?)- 



(A17) 



3=1 



3=1 



APPENDIX B: DISAPPEARANCE OF ODD POWERS OF D IN R 

This appendix proves that the first line of Eq. ( |l9| ) does not generate terms d 2 (dh) 2l+1 , where i is a non- negative 
integer. To this end, consider a Taylor expansion of R l n+k and set the lattice constant to be 0. We call this quantity 
T/ (£)), and the corresponding quantity for S Si (D). One should not confuse this procedure with the continuum 
limit. It is enough to prove that T ; °(— D) = T' ; (_D). It is clear that T Z °(D) has no k dependence because the lattice 
constant is set to 0. This fact enables us to write the recursion relation of T,°(D) in the symmetric form 



2Ti°{D) = SiMTE^D)- 



S lo (D)+D\ fS lo (D)-D 



2 J V 2 
S lo (D)-D\ fS lo (D)+D 



(Bl) 



Hence we see that 



2(T/°(£>) - T/°(-D)) = Si a (D)(T}^ 1 (D) - T^-D)) (B2) 



The mirror symmetry of this system restricts the form of Si to be even in D, so r/ 1 (Z)) = T,^ 1 (— D) implies 

T/°(D) = T/°(-D). Indeed T[°{D) = S h + (Sf o - D 2 )/4 and the logic of the mathematical induction proves the 

disappearance of the odd powers of D in T/°(£>). The vanishing j(l,k) in appendix |X| is a consequence of this 
property. As a result, we can safely affirm that CRSOS model is described by the VLD equation. 
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TABLE I. Numerical values of p(lo), v(lo), and \(lo)- 



lo 


p(lo) 


v(lo) 


A(Jo) 



1 


i 

0.299 027 750 50 


i 

3.3887 x 10" 1 


_ 5 

-9.5709 x 10" 1 


2 


0.293 696 759 81 


9.7310 x 10~ 2 


-2.9204 x 10" 1 


3 


0.292 988 255 45 


1.9569 x 10 -2 


-6.0527 x 10" 2 


4 


0.292 903 676 47 


3.2656 x 10~ 3 


-1.0287 x 10~ 2 


5 


0.292 894 314 95 


4.8299 x 10~ 4 


-1.5406 x 10" 3 


6 


0.292 893 329 94 


6.5674 x 10 -5 


-2.1140 x 10" 4 


7 


0.292 893 229 81 


8.3964 x 10~ 6 


-2.7217 x 10 -5 


8 


0.292 893 219 88 


1.0242 x 10~ 6 


-3.3380 x 10~ 6 


9 


0.292 893 218 92 


1.2038 x 10~ 7 


-3.9408 x 10~ 7 


10 


0.292 893 218 82 


1.3730 x 10~ 8 


-4.5115 x 10" 8 



oo 0.292 893 218 81 



FIG. 1. Relation of species to the height slope. We use the integer to indicate the location of particles and the half-integer 
for the height configuration. Hence the RSOS condition at site n + 1 is determined by the situation at sites n and n + 1. 

FIG. 2. Plots for the saturation width as a function of the system size for various lo- We find the saturation width by a 
least-squares fit and the error bars represent three times the standard deviation, which includes 99% of data. We fit the data 
as a function of L and find two exponents. The exponent for the smaller system sizes are written in the bottom left and that 
of the larger system sizes in the upper right. The lines show the fitting results. Up to 512, it seems plausible to insist that the 
systems are in the scaling regime, but the data for 1024 shows the clear discrepancy within an error bar. 

FIG. 3. Scaling plot of saturated widths of the one-dimensional RSOS/H model. The scaling variable is Iq/L. The values of 
lo are equal to those in Fig. 0. We also draw the line a^w - a kP" = x aA as a guide to the eye. 
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